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The adhesive dynamics of a one-dimensional aggregating gas of point particles is rigorously de- 
scribed. The infinite hierarchy of kinetic equations for the distributions of clusters of nearest neigh- 
bours is shown to be equivalent to a system of two coupled equations for a large class of initial 
conditions. The solution to these nonlinear equations is found by a direct construction of the rele- 
vant probability distributions in the limit of a continuous initial mass distribution. We show that 
those limiting distributions are identical to those of the statistics of shocks in the Burgers turbulence. 
The analysis relies on a mapping on a Brownian motion problem with parabolic constraints. 



I. INTRODUCTION 

Rigorous solutions of the many-body dynamics are extremely rare. The present paper yields up such a precious 
case, reviewing parts of the existing works and completing the research which has been carried out in recent years. 

The system under study is a one-dimensional gas of point particles forming aggregates through perfectly inelastic, 
adhesive collisions. The motion between collisions being free, the process is called ballistic aggregation. This dynam- 
ics is deterministic, randomness occurs only through the distribution of initial datas. In the context of statistical 
mechanics, the model was introduced and studied numerically in Q and further investigated analytically in ||], j^] 
and Q. However in the context of fluid dynamics, it has been recognized much earlier (see e.g. that the 

evolution of shocks in the inviscid limit of the (decaying) one-dimensional Burgers equation obeys the laws of ballis- 
tic aggregation. Moreover when the initial particle velocities are Gaussian and independent (equivalently when the 
initial Burgers held is a white noise) both models are isomorphic to a problem of Brownian motion under parabolic 
constraints. The latter problem has eventually received a solution in closed analytical form Q, Q which enables to 
predict the exact distribution of shocks in the decaying Burgers turbulence, and correspondingly the distribution of 
masses and velocities of the aggregating particles. These close connections between problems of different origins is 
an attractive feature of ballistic aggregation. The purpose of this work is to establish precisely the above mentioned 
relations and to give an exact and fairly complete description of ballistic aggregation. 

There are essentially two approaches. The first one consists in writing the dynamical hierarchy of equations coupling 
C ' the many-particle distribution functions using the tools of kinetic theory, and hopefully solving these equations. The 
second route is by explicitly constructing the dynamics from suitably chosen initial conditions: the construction should 
^ of course yield a solution of the kinetic equations. Both approaches are discussed here and both will make the link 
O with the dynamics of Burgers shock waves manifest. 

O After recalling the way in which one can describe the state of an infinite gas (Section |l|), we present in Section 



[II the derivation of the infinite hierarchy of kinetic equations coupling the time evolution of distributions of many- 



k>( particle clusters of nearest neighbours. The derived hierarchy, after integration over momenta of the aggregates, turns 
' out to be identical to the hierarchy found in the study of the Burgers turbulence |^ . 



In Section IV an original result is derived. The infinite dynamical hierarchy is shown to be compatible with the 
factorization of the many-body distributions into products of two-particle conditional distributions and the one- 
particle density. Such a factorization is propagated by the ballistic aggregation. This yields an exact closure of the 
hierarchy and permits to reduce it to a system of two coupled equations (Section ^) which allow self-similar solutions 
(Section 0). 



*supported by the Swiss National Foundation for Scientific Research. 



1 



The Section VI] is devoted to an explicit construction of the particle distributions. The construction, which has 
already been presented in uses a simple initial condition which enables to reformulate the problem within the 
theory of Brownian motion. To establish the relation with shock wave dynamics at this level, one must envisage the 
delicate question of the continuum limit. Indeed the Burgers velocity field takes values in the continuum, whereas 
masses of aggregates, resulting of sums of elementary initial masses, have discrete values. A strict isomorphism will 
only be obtained in the continuum limit of the aggregation process letting the distribution of initial masses tend to 
a continuous and uniform mass density. In terms of Brownian motion this leads to the subtle problem of controlling 
the cumulated effect of Brownian excursions in infinitely many small time intervals. This problem is addressed and 



extensively discussed in Section VIII : in the continuum limit the particle distributions become identical to those of the 
shocks in Burgers dynamics so that all the results of immediately apply. All distributions of order three and more 



factorize in the way described in Section IV, so giving a complete statistical description of the aggregation process. 
In Section [X the exact predictions for the mass and velocity distributions are compared to existing bounds as well as 
to the findings of previously formulated approximate theories. In particular significant discrepancies with mean field 
theory appear in the domains of large and small masses. The asymptotic form of the two-particle correlation function 
at large distance displays the range of spatial correlations induced by the aggregation process: one finds a very rapid 
decay of such correlations. In the last section, we briefly indicate by which mechanisms our explicit solution verifies 



the closed dynamical equation derived in Section IV (and hence the whole hierarchy), thus completing our analysis 
of the model. 

II. STATES OF THE AGGREGATING GAS 

At any time t > 0, the system is composed of point particles moving in R^. The state of a given particle is entirely 
characterized by specifying its position X, its momentum P, and its mass M. 
Consider a closed space interval [Li, L2] , where Li < £2- We denote by 

fik{l,2,...,k;t\Li,L2) (1) 

the probability density for finding precisely k particles within the interval [Li,L2\ at time i, occupying the sequence 
of states 

j = {X„P,,M,), j = l,2,...,fc (2) 

ordered in space according to the inequalities 

Li < Ai < < . . . < Xfc < L2 (3) 

The set of densities /i/c(l, 2, . . . , A:; t\Li, L2), fc = 0, 1, 2, . . ., defined for arbitrary intervals [Li, L2], provides a complete 
statistical description of the state of the infinite volume of the gas. Summing up the probability weights for all possible 
events within [Li,L2] one gets the normalization condition 



^lo{t\L^:L2)+Y. JdlJd2...Jdk e{Xi-Li)9{L2-Xk) 

k—l 



where 



k 

X n ^(^J - ^l-^)^'k{l, 2, . . . , fc; t\Li,L2) = 1 (4) 
i=2 



+ CXD /' + 00 /' + 00 

dj = I dXj / dPj / dMj , (5) 
) J-00 Jo 



and 9{X) is the unit Heaviside step function. 

Knowing the probability densities 2, . . . , fc; L2) one can evaluate the densities 2, . . . , fc; of the 

nearest neighbours configurations by letting the extremities of the interval [Li, L2] approach the positions of the first 
and the last particle, respectively 

^fc(l,2,...,fc;i) = lim lim ^^(1, 2, . . . , fc; L2) (6) 
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In particular, fii{l;t) represents simply the number density of particles at time t in the one-particle state 1 = 
(Xi,Pi,Mi). 

An alternative way of describing the state of an infinite system consists in determining the complete set of reduced 
distributions 2, . . . , fc; t), representing the number density of ordered k-particle states (1, 2, . . . , fc). The difference 
with respect to the densities /ifc(l,2, . ..,k;t) comes from the fact that the k particles need not represent the set of 
nearest neighbours. Although in evaluating 2, . . . , k;t) the inequalities 

Xi<X2<...<Xk (7) 

are still supposed to hold, one has to consider all configurations of the system compatible with the condition that the 
states (1, 2, . . . , fc) are occupied at time t, with an arbitrary number of particles between pairs {j, j + j = 1, . . . , fc— 1. 
Clearly, there is then no difference at the level of the one-particle densities, and, as it has already been mentioned, 
the equality 

pi(l;t) (8) 

holds. 

For the two-particle reduced density /92(1, 2; i), the corresponding formula reads 

P2(l,2;i) =M2(l,2;t)+ V f dl' f d2' . . . f dr' e{X[ - Xi)9{X2 ~ X',) 



■r-1 



X 



II e{X!^^, - Xj) fi2+r{l, 1', 2', . . . , r', 2; t) (9) 

The possibility of the presence of r particles (r = 0, 1,2, . . .) between the two particles at Xi and at X2, has been 
taken into account in (||). By definition the density /i2+r(l, 1', 2', . . . , r', 2; t) vanishes outside the region Xi < 
X[ < . . . < Xl < X2, reflecting the linear ordering of the particles. 

In a similar way (i.e. by considering all possible intermediate states) one can express pfe(l, 2, . . . , fc; i) in terms of the 
densities of nearest neighbours Hk+ri^, 2, . . . , k + r;t), r — 0,1,2, . . ., for arbitrary fc = 3, 4, . . .. The general relations 
between both types of description of the states of an infinite system are discussed in Ruelle's book [bI. 



III. EVOLUTION OF THE DISTRIBUTION OF NEAREST NEIGHBOURS 



Our main object here is the study of the dynamics of aggregation. The point masses forming the system move 
freely between collisions. When two of them collide, they instantaneously merge forming a new particle whose mass is 
equal to the sum of their masses. As the momentum is also conserved, the aggregate continues the motion along the 
center of mass trajectory of the pair it has been formed from. Hence, the microscopic dynamics consists of periods of 
a free motion separated by perfectly inelastic (sticky or adhesive) binary collisions. 

Suppose that at some moment two particles occupy the states a = {Xa, Pa, Ma) and b = (Xb, Pf,, Aft,), with Xa < Xf,. 
The notation (a -I- b) will be used to denote the state corresponding to the center of mass motion of such an ordered 
pair. Hence 

We shall restrict our analysis to spatially homogeneous systems. The one-particle density does not depend in this 
case on the position variable 

/ii(l;t) = ^i(Pi,A/i;i) (11) 

Therefore, the changes of /Lti(l;i) in the course of time are exclusively due to collisions. Clearly, in one dimension 
only the nearest neighbours can collide. The rate of collisions between adjacent particles j and {j + 1) is proportional 
to the density of the nearest neighbours fj.2ij,j + 1;^)- It can be written as 

C{j,j + i)fi2{j,j + i;t) (12) 

where the collision factor 
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C(j,j + 1) 



6iX,+i - X, - 0- 



(13) 



chooses (with the help of the Dirac (5-distribution) only the precolhsional configurations. As the motion between 
colhsions is free, the rate ( [l^ ) is proportional to the relative velocity of the colliding pair (the 0-factor in C{j,j + 1) 
assures the mutual approach of the particles). 

In order to evaluate the time derivative of the one-particle density ( [Til ) let us consider first the events which lead 
to the creation of the state 1 by adhesive collisions. Their rate of occurrence is given by 



dl' / d2' C{1', 2')S[1 - (1' + 2')]fi{l', 2'; t) 



Here the distribution 

S[l - {!' + 2')] = 5 



Xi- 



M[X[ + M^X^ 



5[P, - {Pi + P^)]S[Mi - {M[ + M^)] 



(14) 



(15) 



picks out only those collisions which create the aggregate in the state 1. 

Consider now the events leading to the destruction of state 1. When the precolhsional state of one of the particles 
is 1, the collision removes this state from the system. The rate of such annihilating events is given by 

(16) 



di'C(i',i)A*2(i',i;t) + j di'C(i,i')M2(i,i';i) 

Combining the gain term ( p^ ) and the loss terms (p^), we arrive at the equation 
-Hi{l-,t) = j dl' j d2'C(l',2')<5[l-(l' + 2')]/i2(l',2';t) 



d_ 
di' 



- / dl'C{l'^)^l2{l'^\t)- / di'C(i,i')A^2(i,i';i) 



(17) 



In a similar way one can derive the evolution equation for the two-particle distribution of the nearest neighbours. 
It reads 



^+Li2 + C(l,2) 



M2(l,2;t) = 

dV / d2'C{l\ 2'){5[1 - {!' + 2')]a^3(1', 2', 2; t) + 5[2 - {!' + 2')]a^3(1, 1', 2'; t)} 
dl'C7(l',l)/i3(l',l,2;t)- / dl'C7(2,l')Ai3(l,2,l';i) 



where 



-^12 — 



Pi d 



P2 d 



Ml dXi M2 dX2 



(18) 



(19) 



is the generator of free streaming. On the left hand side there appears the loss term C(l, 2)/i2(l, 2; t) describing the 
destroying effect of a possible merging of the pair (1,2). The right hand side takes into account the processes of 
creation of aggregates in the states 1 or 2 (the first two terms). Finally, the last two terms represent the destruction 
of the two-particle state (1,2) by collisions with the left nearest neighbour of particle 1, and with the right nearest 
neighbour of particle 2, respectively. Equation ( p^ ) relates the rate of change of /i2(l, 2; i) to the three-particle density 
Ai3(l,2,3;t). 

Denoting by ii...A; the k-particle generator of free streaming 



Mi dXi 



(20) 



we write the general /c-th equation of the infinite hierarchy (/c = 1, 2, . . .) in the form 
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fc-1 



Aife(l,2,...,fc;t) 



iZh^' I d2'C(l',2')<5[j-(l' + 2')]Mfe+i(l,...,J-l,l',2',j + l,...,fc;i) 
- j dl'C(l',lK+i(l',l,2,...,A:;t)- ^ dl'C(fc, 2, ■ • ■ , A:, 1'; t). 



(21) 



The fact that /ifc(l, . . . , j, j ' + 1, . . . , fc; is the density of a sequence of nearest neighbours imposes a restrictive 
condition on its support. Indeed, tracing backward in time the free trajectory of particle j to the moment t = 0, we 
find the position of the center of mass of the part of the initial system it has been formed from. So, if we denote by 
p the constant mass density at t = 0, the point 



X 



right 



{-t) = X, - 



El 



t + 



M. 

Yp 



3 — 



(22) 



has the meaning of the right extremity of the initial mass region which produced the mass Mj through adhesive 
collisions. This right extremity must coincide with the left extremity 



p. 



M. 



M. 



2p 



X, 



(23) 



of the region which produced the mass A/j+i. Otherwise, the pair (j, j + 1) would not represent nearest neighbours. 
In (|2^) and ( p3| ) we have introduced the notation 



El. 

Mi 



-t± 



2p 



We thus conclude that the density /i/c(l, . . . ,k;t) contains necessarily the singular factor 



fc-i 



l[S[X,+,-X,~Y+,+Yf] 



(24) 



(25) 



This conclusion implies important consequences for coUisional configurations where Xj+i = Xj + 0+. The distribution 
(p5|) imposes then the condition 



^+1 = 

Hence, in the hierarchy equations ( pT| ) one can replace the collision terms C{j,j + 1) (see (|l3|)) by 



2pt 



-5{X,+, - X, - 0- 



(26) 



(27) 



With the use of (|27|), one finds, after integrating over all momenta, the infinite hierarchy for the distribution of masses 
and distances between the particles which has already been derived in the study of the Burgers model of turbulence 
1^ . Our independent reasoning based on the microscopic laws of the aggregation process confirms the existence of a 
one-to-one correspondence between the two problems. 

It is quite remarkable that the infinite set of equations ( ^l| ) can be rigorously reduced to a system of two coupled 
nonlinear equations under a simple assumption, which can be shown to be satisfied for a large class of initial conditions. 
This important fact, derived in the next section, has not been noticed in The authors restricted their analysis 
therein to the first equation of the hierarchy, supplemented with an integral relation obtained by integrating equation 
( p^ ) over both momenta and masses. In the search for self-similar solutions, they assumed a specific form of two-particle 
correlations, which eventually led them to an erroneous conclusion that the masses were exponentially distributed 
. The exponential distribution has been eighteen years later conjectured again in and shown to follow from 
the hierarchy equations (|2^) within the weak mean-field approximation j^]. It turned out that a rigorous analysis 
invalidated these approximate results (see 0] and section IX). 
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IV. REDUCTION OF THE HIERARCHY 



A fundamental role in further considerations will be played by the conditional probability density 

M2|i;0 = ^^^%^ (28) 

for finding the right nearest neighbour of a particle supposed to occupy the state 1, in the state 2 at time t. The 
normalization 



J d2 Ai(2|l;i) = 1 (29) 



expresses the fact that in a homogeneous system the right nearest neighbour does exist in some state with certainty. 
The singular factor 

5 [X2 - - Y+ + Ff ] (30) 

present in ^2(1, 2; t) (see (|5|)), and thus also in /z(2|l; t), involves the combination Y{~ = tPi/Mi~ Mi/2p of mass and 
momentum of particle 1. This combination appears in formula ( |2^ ) for the position of point X"^^^ {—t) , separating 
from the right hand side the masses which contributed to the formation of particle j from the rest of the system. 

The reduction of the dynamical hierarchy ( pT| ) presented in this section is possible if /i(2|l;t) depends on Pi and 
Ml exclusively via the effective variable . The following structure of the conditional probability density will be 
thus assumed 

Ai(2|l; t) =d[X2-Xi- Y+ + Ff ] fl{M2, P2, i^f ; (31) 



In fact, it will be shown by construction in Section VII that for a large class of initial conditions the above hypothesis 
is verified. 

The presence of factor ( p5| ) in /^^(l, . . . , k; t) clearly indicates the existence of correlations between the states of the 
nearest neighbours. Let us suppose that only two-particle correlations exist, so that 

j=k-i 

/Xfc(l,2,...,fc;i) =Mi(l;i) n /^(J' + Ib';^) (32) 

i=i 

We shall prove now that the above factorization is compatible with the dynamics of the system, and that the factorized 
densities (^) yield a solution to the hierarchy equations (^l|), provided ^,i{l\t) and /i(2|l;t) satisfy the system of 
coupled equations 

^/.i(l;t) = j dV j d2'C(l',2')<5[l-(l' + 2>i(lV)M2'|l';i) 

dl'C{l\l)ni{l'-t)n{l\l'-t)- f dl'Cil,l')fii{l;t)fi{l'\l;t) (33) 



-+Li2 + C(l,2) 



/i(2|l;i) = J dV J d2'C(l',2')5[2-(l' + 2')Ml'|l;i)M2'|l';0 
dl'C(l,l')M(l'|l;t)- / rfl'C(2,l')A^(l'|2;<)lM(2|l;t) 



(34) 



Equation ( ^3[ ) is just the first equation of the hierarchy ( p7| ) in which /i2(l, 2; t) has been replaced by /ii(l; i)/i(2|l; t). 

Let us consider the general equation ( pi] ) for k > 2. Upon inserting the factorized form (^) of the densities /i^ and 
l^k+i we get (with the use of the first equation (pS)) the following relation 



Mi(i;0 



d_ 
dt 



fc-i 

ii...fc + 5]C(j, j + 1) 



fe-i 



\{^l{r + l\r■,t) 



r=l 
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/ dl' [ d2'C(l',2')<5[l-(l' + 2')]Ai2(l',2';t)[M2|2';t)-M2|l;t)] \{ + 

k J ~2 

/ dV j d2'C{l\2')5[j ~ {!' + 2')]\{^,{r +l\r-t) 



k-1 



x^x{l'\]-l■t)^l{2'\l'■t)^Ji{J + l\2'■t) \[ ^Ji{s + l\s■t) 

k-1 



-^l^{l■,t) 



di'c(i,i')/i(i'|i;0 " / d\'c{k,l')^JL{l'\k■t) 



The first term on the right hand side vanishes. Indeed, equation (^6|) imphes the equaHty 

Ml 2p M{ + 2p 2p 



(35) 



(36) 



So, taking into account the supposed structure ( ^l|) of conditional densities, we can replace /i(2|2'; t) in (|35| ) by 
^(2|l;t). Using exactly the same reasoning one can show that in the remaining gain terms in (pq), involving the 
creation of particles j = 2, 3, . . . , A: through aggregation of 1' and 2', the density /i(j + 1|2'; t) equals p,{i + t). This 
permits to rewrite the hierarchy equation (pq) in the following form 



d_ 
di 



k-l 



ii...fc + ^C(j,j + l) 



k J ~ 2 k—\ 

Y.jdi'j d2'c(i',2')%-(i' + 2')]nM(r + iht)M(i'b--i;%(2'|i';t)nM^^ + i|^;^) 



dl'C(l,lXl'|l;i) -y (il'C(/e,l')Ai(l'|A:;i) 
The last term in (B7|) can be conveniently rewritten as 



fc-i 



(37) 



dl'C{j~l,l')p{l'\]~l-t)- I dl'C{3^')^l{l'\3■t) 



E 

i=2 



[]Ai(r + l|r;t) 



(38) 



A straightforward calculation shows then that equation (^5|) is satisfied (for any /c > 1) if the conditional probability 
density /i is a solution of equation (^4|). 

We have thus proved that solving the system of coupled equations ( ^3| ) and (^), one obtains the solution of the 
infinite hierarchy (|2|) in the factorized form (^). This remarkable reduction of the dynamical hierarchy to a closed 
system of two equations does not involve any approximation, provided ^(2|1; t) has the structure (^). The adhesive 
collisions do correlate the motion of the nearest neighbours (the mean field approach is ruled out), but no three- or 
more-particle correlations are created in the course of time. The same kind of situation has been already discovered 
in the case of ballistic annihilation |10| . 



V. EVOLUTION OF REDUCED DENSITIES 

The factorized form (^2|) of distributions pk implies the analogous factorization of the reduced densities 

p(i,...,fc;i) -Pi(i;t) n -"(J + Ili;*) 



(39) 



where p(2|l;t) is the conditional density 
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p(2|l;t) EE p(X2i,P2,M2|0,Pi,A/i;0 = ^^^^ (40) 



From the closed system of equations (^3|) , (^J) for the nearest neighbour distributions one can derive the corresponding 
set of dynamical equations for the densities pi = and p(2|l; t). First of all let us notice that the relation implies 
the equality 

P2 {Xi , Pi , Ah , P2 , Afa; t) = ^^2 (^1 , Pi, Ah, Xi+, P2, Ah; t) (41) 
Hence, equation ( p3| ) preserves its form 

|pi(l;i)= / dl' /''^2'C(l',2')5[l-(l' + 2')]pi(l';i)p(2'|l';t) 



di'c(i',i)pi(i';Op(i|i';i)-y rfi'c(i,i')pi(i;Op(i'|i;0 (42) 

In order to derive the closed evolution equation for p{2\l;t) one has to use the formula 
p(2|l; i) = fi{2\l; i) + X] / / ' ' ' / ^^^^ " ^i)^^(^2 - ^;) 



X 



n 0(X;.+i - X;) M(l'|l; t)p{2'\l'; t) . . . /.(r'l/ - 1; t)p{2\r'; t) (43) 



(see (^ and (^)). Using then repeatedly equation (^) one arrives at the equation 
^ + Li2 + C(l,2) 



p(2|l;<) = J dV j d2'C(l',2')<5[2-(l' + 2')Hl'|l;i)p(2'|l';t) 
dl'C(l,l')p(l'|l;t)b(2|l;t)-p(2|l';t)] 

dl'{C(l',2)p(l'|l;t)p(2|l';t) + C(2,l')p(l'|2;i)p(2|l;0} (44) 



Equations (0) and ( ^4[ ) form a closed system which suffices to determine the reduced densities of any order owing to 
the relation ((391). 

VI. SELF-SIMILAR SOLUTIONS 

It seems interesting to check whether the aggregation dynamics is compatible with self-similar solutions. In other 
words, whether the evolution of the merging masses can be entirely reduced to rescaling of densities of a given shape 
(this idea has been followed in the study of the Burgers model of turbulence [^). In order to investigate this question 
we insert into equations (|3^), (^ ) the assumed scaling formulas 

Hi{Pi,Ah;t) = pi {Pit^-^'',Aht-''; 1) (45) 
/x(X2i, P2, M2IO, Pi,Ah;t) = t-<' ^(X2it^", P2t^-2", M2t""|0, Piti-2", Afii-"; 1) (46) 

where X21 = X2 — Xi and where we have used the fact that the three variables 

X21, —t, — (47) 
AI ' 2p ^ ' 

have the same dimension and must thus scale in the same way. 
As the mass is conserved, the integral 



/+00 />oo 
dP dAdAdf^' i^{Pt^~^",Me) 
-00 Jo 



(48) 
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representing the total mass density does not depend on time. This implies the relation 



71 = 1 - 4a. (49) 
On the other hand, the normalization condition ( |2g| ) leads to 

72 = 1 - 4a. (50) 

It can be checked by a straightforward calculation that the distributions /ii(lji) and /i(2|l; t) of the form ( ^ ) and 
(^), respectively, with the exponents (a, 71, 72) satisfying equations (^ ) and (po|), lead to a consistent closed system 
of equations for the self-similar distributions /xi(P',M'; 1) and ^(2'|1'; 1) when put into the reduced hierarchy (|3^), 
(p^). It follows that the dynamics does not suffice to fix completely the values of the exponents. They can thus 
depend on the nature of the initial condition of the system. 

VII. A STATISTICAL MECHANICAL MODEL 

In this section, we give a microscopic description of the dynamics of ballistic aggregation from the view-point of 
statistical mechanics by recalling the model introduced in and |^ . 

One considers the aggregation process developing from a simple initial condition. At t — 0, all the particles 
have the same mass m and are located on the sites Xj = ja of an infinite regular lattice with lattice constant a 
(j = 0, ±1, ±2, ...). The mass density p — m/a. It is assumed that the initial momenta are uncorrelated, with a 
distribution corresponding to thermal equilibrium at inverse temperature /3: 



As it has already been explained in Section III, the dynamics of the sticky gas has the remarkable property that any 
mass aggregate is found on the trajectory of the center of mass of the initial cluster it has been formed from. In 
particular, the state of an aggregate Mj — n.jm at time < > determines uniquely the set of Uj consecutive initial 
masses which constitute it: they are located at i = within the interval \Xj — -I- a/2, — Y~ — a/2] (see also 
Eqs.(p2|),(l2^)). This property permits to write explicitly the kinematical constraints that select the subset of initial 
phase space configurations leading to the occurrence of a specific sequence of aggregate states j = {Xj,Pj,Mj), j = 
1, . . . , A:, within [Li, L2] at time t > 0. The distributions /ifc(l, 2, . . . , /c; t \Li, L2) are then obtained by averaging these 
constraints over the initial state. 

The origin of the constraints is two- fold. One class of them ensures that the rij initial particles located in the interval 
[Xj—Y^ + a/2, Xj — Y~ — a/2\ do merge and form the aggregate before time t. A necessary and sufficient condition 
here is that the trajectories of the centers of mass of all the pairs of subclusters of consecutive initial particles, which 
form a partition of the initial rij -particle cluster, cross before time t. After averaging, this leads to functions denoted 
by ImiMj, Pj]t) which are the probability densities for the formation of masses Mj = njin, with momenta Pj from 
the corresponding sets of Uj initial neighbouring masses. 

The second class of constraints guarantees that no particles other than the k specified aggregates are found within 
[Li,L2] at time t. The particles initially located in {— 00, ja] will be found to the left of Li if the positions of all 
the centers of mass of clusters of consecutive initial masses located in [{j — i)a,ja\, i = 1, 2, . . . stay smaller than Li 
up to time t. After averaging, this will be expressed by a function Jjn(Y;t), where J„i{—Li + (j + l/2)a;t) (resp. 
Jm{L2 — {j + l/2)a;t) is the probability for finding the particles initially located in (— cx3, ja] (resp. [(j -I- l)a,oo)) in 
(— cx),Li] (resp. in [L2,oo)). In particular, in terms of function defined in this way the probability /^o,m(^|-^i, ^2) 
of finding the interval [Li, L2] void of particles is given by 

/^o,m(i|ii, ^2) = ^ Jm i-Li + {j + l/2)a; t) J,„ (L2 - (j + l/2)a; t) (52) 
j 

The probability density fii^m (^ reads |^ 



'^We add the index m to keep in mind the discrete nature of the initial state. 
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Ml,' 



J2 SiMi~mm)Y,S{X,-Y+ -{j + l/2)c 



ni — 1 



X {Ah ,Pi;t) J™ {~Y+; t) J„ (Ff ; (53) 



where the summation involving ^-functions is a manifestation of the discreteness of the initial masses and the lattice 
positions. The function in ( p3| ) ensures that the particles initially located in [Xi — + a/2, Xi — Y{~ — a/2] have 
met to form the aggregate, while the functions Jm ensure that all the particles initially located on the left of Xi — Y^ 
(resp. on the right of Xi — Y^) are at time t on the left (resp. on the right) of Xi. 
The probability density /i2,m (^ is given by 

oo oo 

M2,™(l,2;^) = ^?(X2-Xl) J2 5(Mi-n^m) ^ SiM^-n^m) 

rii— 1 ri2 — 1 

X ^ (5 (Xi - Y+ - (j + l/2)a) 5 {X^ - X^ + Ff - Y+) 

j 

X /„ (Ml , Pi ; t)Ira (M2 , P2 ; t) Jm {-Y+ ; t) Jm [Y^ t) . (54) 

In this expression the jii and j summations involving (5-functions reflect again the discreteness of the initial conditions. 
The additional (5— function obtained here by construction reflects the fact that the aggregates are nearest neighbours 
and imposes precisely the condition ( ^l| ) under which the hierarchy could be reduced. Once again, functions Im 
stand for the formation of aggregates 1 and 2, while functions Jm ensure that particles initially located on the left of 
Xi — Y^ (resp. on the right of X2 — 1^2") ^^^y time t on the left of Xi (resp. on the right of X2). 

The higher order distributions ixk.m can be constructed along the same lines. They verify the factorization (p2) 



/Xfc,m(l, . . . ,fc;t) = Aii,m(l;i) Y\ Mm(j + (55) 



with the conditional probability /im(2|l;i), derived from ( p3| ) and (|5^), of the form 

^2,m(l,2;i) 



/^m(2|l;i) 



Aii,m(l;i) 



°° J iY~-t] 

= e{X2 - Xi) ^ 8{M2 - n2m)5 {X2 - Xi + Ff - Y+) ImiAh, P2]t) /" ' ' (56) 

712=1 

Note that this conditional probability has the structure (pl|). 

To have the model in an explicit form, it remains to give the formulae that express functions Im and Jm in 
terms of the constraints: the details of the calculation can be found in The important point is that, due to the 
uncorrelated Gaussian initial velocity distribution, the whole aggregation dynamics can be mapped on the following 
equivalent Brownian motion problem: 

Let P{t) be a Brownian path starting from the origin, P(0) = 0. Then, for M — nm, 

Im{M,P;t) = E(^n^o){P{rm) > rm[{M ~ rm)/2pt + P/M],r ^ 1, . . . ,n- l|P(nm) P} (57) 

is the conditional measure of such paths constrained to be above the parabolic barrier rniKM — rm)/2pt + P/M] at 
discrete "times" Tr = rm, r = 1, . . . ,n — 1 and to end in P at "time" nm, see Fig.|[ 
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FIG. 1. Brownian interpretation of the function Im{M, P;t). The Brownian motion starts at P(0) = 0, ends at P{M) — P 
while overpassing the points P{rm) > rm{2pQ^ — rm)/2pt, (r = 1, . . . ,n — 1) with Q"*" = Pt/M + M/2p. 

Likewise 

J^{Y- 1) = £;(o,o){P(rm) > rm{2pY - rm)l2pt, r = 1, 2, . . .} (58) 

is the measure of the paths that remain above the barrier rfn{Y — rm)/2pt for all discrete times = rm, r — I, . . 
see Fig.|^. We recall that the dynamics is deterministic, and randomness enters only through the initial velocity 
distribution: P(t) is a process in momentum space, as a function of the mass of aggregates. 




FIG. 2. Brownian interpretation of the function Jm{Y;t). The Brownian motion starts at P(0) = and overpasses the 
points P{rm) > rm{2pY — rm)/2pt, {r — 1,2, . . .). 

A simple and immediate consequence are the scaling properties of functions and Jm with respect to time. Owing 
to the fact that the scaled Brownian motion {pt)'^/'^P{T{pt)~") is equivalent in probability to P{t), the functions 
and Jra obey the scaling relations with a — 2/?, 

I^{M,P-t) = (pt)-i/3j^„,(Af',P';l) = {pt)-^'''lra'{M',P') 

JmiY;t)^J^.iY';l) = Jm'{Y') (59) 

with 

m' = m(p^)-^/^ P' = P{pty^'^, M' = nm' = M{pt)''^/^, Y' = {pty^^^pY (60) 

The construction (p3[)-(p8|) yields in principle a solution to the dynamical hierarchy that has the factorization property: 
Mi,m(l;i) and /im(2|l;t) have to verify the coupled equations (p3|)-(p4), adapted to the case of discrete masses. We 
shall not proceed to this verification now, but rather simplify first the discussion by taking the continuum limit of our 
statistical mechanical model. 
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The continuum Hmit amounts to let to ^ and a —^ while keeping the initial mass density p = m/a and time 
fixed. Because of (^9|) one can alternatively fix the scaled momentum P' and the scaled mass M' and look for the 
large time asymptotics t — > oo. In the first view, one looks for the formation of masses of order 1 at a given time 
arising from an initial infinitesimal dust. In the second view one looks for the formation of masses of size ^ mt^^^, 
t —^ oo, while keeping the discrete initial condition. These two views are equivalent. The interesting point is that 
in the continuum limit the Brownian expressions ( ^7|) and ( ^8|) can be computed in a closed analytical form, thus 
providing an explicit complete statistical description of the aggregation process. Moreover this description exactly 
coincides with that of the statistics of shocks in the inviscid limit of the Burgers equation with white noise initial 
data. 

To conclude one can check that the factorization property (^5|) persists for a wider class of initial conditions, 
for instance allowing any uncorrelated initial velocity distribution, and an arbitrary (non random) choice of initial 
positions and masses. 



VIII. THE CONTINUUM LIMIT 



In this section, we determine the continuum limits (to — > with p = m/a fixed) of the one- and two-point probability 
densities 

lim pi^m{^;t) = /ii(l;t) 

m— >0 

lim ^2,m(l,2;t) -A*2(l,2;t). (61) 

m — »U 

Owing to the scaling relations ( |59| ) it is sufficient to calculate functions pi{l;t) and /i2(l,2;f) for t = 1 andQ P = 1- 
From now on we drop the time parameter from the notation setting simply /ii.m(l;t) = /xi.m(l), /i2,m(l,2;t) — 
/i2,m(l, 2), and so on. Later in this section we will comment on the way to recover the time variable in the contimmm 
limit of these densities. 

Let us first introduce the transition probability density kernel 

KmAMi,Pi,M2, P2) = E^M,,P,){P{rm) > Mrm), r = m + 1, . . . , - l|P(Af2) = P2} (62) 

for the Brownian motion P(r) to start at P(Mi) — Pi, and end at P{M2) — P2 with Mi — riim and M2 = n2m, 
while overpassing the points f^{rm), r = ni + 1, . . . , 712 — 1, where fuiT) is the parabola 

Mr)^'^r-^. (63) 

According to definitions ( p7| ) and ( ^8| ) , the functions Im and J™ appearing in the one- and two-point densities can be 
expressed in terms of the transition kernel (p33) as 



I„,{M, P) = K^^Y^ (0, 0, M, P) (64) 

and 

/•oo 

J„(r)= lim / dPi^™.y(0,0,M,P). (65) 

In (|62|), Brownian paths are allowed to make excursions through holes of width to separating the discrete points 
rm. In the continuum limit to — > 0, the weight of such excursions becomes vanishingly small and the paths become 
constrained to overpass the continuous barrier fuir), Mi < t < A/2. Thus for Pi > f,j{Mi) and P2 > fiy{M2) 

K^AMi,PuM2, P2) = /i,(Mi, Pi, M2, P2) + P™,.(Mi, Pi, M2, P2) (66) 

where limm^o Rm,v{Mi, Pi, M2, P2) =0, and 

K,{Mi,Pi,M2, P2) - P(Mi.Pi){^M > fAT),Mi <T< M2\P{M2) - P2} (67) 



^In Q the density was set equal to 1/2. 
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is the transition kernel for a Brownian motion with a continuous parabohc barrier. It satisfies the diffusion equation 



dM- 



2(3 9P5 



■K,{Mi,Pi,M2,P2 



(68) 



with K^{M,Pi,M,P2) = 5{Pi - P2), and the Dirichlet conditions on the barrier: i^^(Mi, Pi, M2, P2) = when 
Pi = MMi) or P2 = /.(M2). 

Our object now is to express the distributions of aggregates in the continuum limit in terms of the kernel (^7|), which 
can be explicitly computed (Section II of |^|). However, in the continuum limit of /,„ and J„i the starting point (0, 0) 
of the Brownian motion lies on the parabola ( |6^ ) where K^{0, 0, M, P) — 0. Thus the determination of the densities 
/ii(l) and /i2(l,2) requires the evaluation of the leading term in the asymptotic expansion of Km,^{0,0, M, P) as 
TO ^ 0. This is not an elementary task since it involves the control of the cumulated effect of Brownian excursions in 
small intervals ((r — 1)to, tto), r = 1, 2, ... in the neighbourhood of the origin. The result is given in the following 
proposition 



Proposition 



lim 

m — >0 



if„^,(0,0,M, P) 



to 



^4i^.(Mo,Po,M,P) 



M > 0, P > U{M) 



(69) 



(Mo,Po) = (0,0) 



With this result the continuum limits to 
initial mass density p — m/a, read 



lim 

m^O 



of functions /„i and Jm, as expressed in (|64|) and (|6^), with fixed 

/„,(M,P) _ 1 



2/3 dPidP2 



and 



lim 



lim 



Ky+{0,Pi,M,P2] 



d 



I{M,P) 



Pl=0,P2 = P 



dP —Ky{Q,PuM,P) 



JiY). 



(70) 



(71) 



Pi=0 



In obtaining ( |70| ) we have taken into account that in (|6^), in addition to (0, 0), the paths have a second contact point 
with the parabola at (M, P) . 

Then, using equations (|53|) , ( 54 ) , the continuum limits of the one- and two-point probability densities can be readily 
found. One multiplies and divides (^^ by and takes (|7^) and (|7^) into account. Then, to ^ playing the role of 
an infinitesimal, the discrete sums go to the corresponding integrals, which can be evaluated owing to the (5-functions, 
yielding eventually the formula 



Mi(i) = /(Afi,Pi)j(-y+) j(rf) 

and, in the same way 

/X2(l, 2) = 9{X2 ~ Xi)S {X2 -Xi+ Yf - Y+) /(Pi, Mi)/(P2, M2) J (-^1+) J {Y^' 



(72) 



(73) 



The time and the initial density dependence can always be reintroduced in the continuum limit. Indeed, through 
Eqs.(M) and (uG),(fl]) we get the scaling properties of functions / and J 



J{Y; t) = (pt)-i/3 j(r'), /(M, P; t) = {pt)-'l{M\ P') 



(74) 



with the scaling functions J{Y) and I{M, P) defined above. This implies the scaling behavior of the densities in the 
continuum limit of the form 



lim pk,rn{l, . . . , /c; t) = /Xfc(l, . . . , fc; t) = /(pt)-5'=/Vfc(l', ■ ■ ■ ,k') 

m— i-0 



with 



/ = (X',P',AP') 



/ pX, 



P. 



(75) 



(76) 
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The scahng functions /ii(l) and /i2(l,2) are given by ( [72| ) and ( [73[ ) and, for higher orders, the scahng functions are 
given by 



fe-i 



(l,...,fc)=/ii(l)n 



A^2(i,J + 1) 
Aii(j) 



(77) 



From Eq.(|52|), wc find that the probabihty density to find an interval [0, x) void of aggregates scales in the continuum 
limit as 



with x' = px{pt) ^/"^ and 



fJ.o{x;t) = ^iq{x') 



fioix) = / dyJ{y)J{x - y). 



(78) 



(79) 



Since the functions I{M, P) and JiY) have been explicitly computed in ||^, the results ( [72|j7^ ) and ( [TT] ) give a complete 
solution of the ballistic aggregation model in a closed analytical form. We pos tpone the discussion of this solution to 
the next section and devote the rest of the present section to the proof of (|6^) . 



Proof of the proposition 

The strategy to prove Eq.(|6^) is to bound the transition kernel ( |62[ ) while linearizing the constraints imposed on 
the first k = M^/m points, < Mq < M. 
We have 



y-'^]T< Mt) <iyT, < r < Mo 



Let us introduce the kernel 

Kl^^,{0, 0, Afo, Po) - ^(o,o){^(™i) > '^rm, r = 1, . . . , fc - l|F(Afo) = ^o} 



(80) 



(81) 



for paths starting at (0,0), ending at (Mn,Po)i while overcoming the linearly distributed sequence of discrete points 
i>rm, r = 1, . . . , A: — 1. Clearly (see Fig. 



Kl^^,{0,0,Mo,Po) < K^AO,0,Mo,Po) < i^,'„^,„Mo/2(0, 0, Afo, Po)- 



(82) 



P (T) 




FIG. 3. Illustration of the "linearized" bounds used to in (B3). 
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Using the Markov property of the Brownian motion, one infers from (^2|) the bounds for the kernel JC,„^,y(0, 0, M, P) 
for M > Mo > 



with 



and 



i^(<) (0, 0, M, P) < KraA^, 0, M, P) < K^liO, 0, Af, P) 
(0, 0, M, P)^ J dPo Kl^^O, 0, Afo, Po)K^AMo, Po, A/, 

K^IUO, 0, A/, P) = / dPo Al,,_Mo/2(0, 0, Afo, Po)ifm,.(Mo, ^o, M, P). 



(83) 



(84) 



(85) 



We shall calculate these bounds by first letting to ^ and then Mq 0. They will be shown to coincide in this limit. 

It is convenient to relate the transition kernel K^^ ^, (^l|) for the Brownian process with discrete linear barrier 
P(rm) > vrm, r — 1, . . . , A: — 1 to that of a dimcnsionless Brownian process ^(t) with covariance equal to 1, the paths 
constrained to be positive at integer times, and with the corresponding transition kernel 



Gk{q) = i^(o,o) {q[n) > 0, n = 1, . . . /c - l\q{k) = q} 



(86) 



For this later process, we define the average {f{q))^. = J dqGk{q)f{q)- The moments {q-')^. , j = 0, 1, 2, . . . will be for 
us of particular interest. 

The relation between the two processes is given by removing a linear drift 



which leads to the relation 



/r^„.,(0,0,Mo,Po) = exp 





{P{t) - vt) 



(87) 




(Po - i^Mo) 



When this is inserted into (|84|), one obtains 



Ki<l{0,0,M,P)^ 



exp 



X Jdq Gfc(q)e-'^(2m;3)V2,^^^^ j^^^^^ ^^^^ ^ ^J2m^^ 



(89) 



Then, we take the following steps. First, for A/q > we can replace in ( p9D the kernel K,n,v for the discrete parabolic 
barrier by its continuous limit (|6^). Next, we expand the integrand up to second order in the variable \/rnq. This 
gives 



i^(<)(0,0,A./,P)=e-'^''^*^o/2 



(l),X„(A'/o,i^A/o,A.f,P) 



+^^^ (g>fc ( -y^K,{Mo, i^Mo, Af, P) + 



+ (<z2)^0(m) + (l),o(TO). 



2 d 



K,{Mo,Po,M, P)|p„ 



(90) 



Recalling that k — Afp/w, we have to evaluate the first moments ('7-')^, of the distribution Gk{q) for large k. This is 
done in appendix ^ with the help of the Sparre- Andersen theorem together with a Tauberian theorem. One gets 



(91) 



, , Mr, , , 1 , 9> I TO 

^ V ttto' ^^'^ 2' V TrAfo 



Hence 
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lim 

m^O 



K,{Mo,vMo,M,P) 



+ - -iy^2f3K,{Mo, lyMo, M, P) + 



2 d 



K4Mo,Po,M, P)\Po=.M, 



Finally, since K^{Mo, Pa, M, P) is differentiable with respect to Pq and vanishes at Pq 
K^{Mo, vMo, M, P) = ©(Mo), A/q ^ 0, so that 



lim lim 



i^i<i(0,0,M, P) 



(92) 

0, one has 
(93) 



(Mo,Po) = (0,0) 



The upper bound Km]^{0, 0, M, P) is treated in the same way and attains the same limit, thus leading to the result 
Eq.(|69[) of the proposition. 

The analysis performed in B was based on the lower bound 



K^.,{0,0,M,P) ^ 1 



m 



d 



27r/3 dPo 



K,{Mo,Po,M, P) 



(94) 



(Mo,Po) = (0,0) 



obtained by retaining the effects of Brownian excursions only in the first interval (0, to) at the origin. The result 
obtained therein differs from the exact limit ( |69| ) by a factor I/a/tt- In fact, the limit ( |69|) involves contributions of 
infinitely many intervals, which can be summed up by applying the Sparre- Andersen theorem. 



IX. EXPLICIT SOLUTION 

The function I{M,P) ([70|), explicitly computed in j^], has the form 



where we set b = (/3/2)^/'^ and 



/(Af, P) = 2b^ exp -6-^ 



X(A/) = ^e- 
fc>i 





■p2 


Af3- 






M " 




) 



X(Af), 



The function J{Y) {\m reads 



with 



J{Y) = N/te-»"'""/»j(y) 



dw 

ioo Al(w7) 



In (|9^) and (psD, Ai(?«) is the Airy function solving the differential equation 

S"{w) - wf{w) = 0, 



(95) 



(96) 



(97) 



(98) 



(99) 



Ai(iy) is analytic in the complex w plane, and has an infinite countable set of zeros — Wfc on the negative real axis, 
< wi < tt'2 < • ■ ■• The asymptotic behavior of these functions has been derived in m 



J(/i) 



, /i — !■ 0; I(/i) ~ exp(— a;i6/i), fj, ^ oo 



and 



-buLJi 



J{U) 



Ai'(-a;i)' 



u — !■ oo] J^{u) ^ —2buexp 



(100) 



(101) 
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In fact, these results, which have been announced in Q, were derived in details in the related framework of the 
Burgers turbulence . The solutions of the one-dimensional Burgers equation with a white-noise initial condition, 
develop, in the limit of vanishing viscosity, a train of shock waves. In Burgers theory, a shock located at X is 
characterized by two parameters /i and 77, and can be identified with a particle of mass M — fi at point X with 
momentum ^ P = —rj. Then it turns out that dynamics of shocks is completely equivalent to ballistic aggregation 
subject to mass and momentum conservation Moreover, the white noise covariance is D/2 = 1/(2/3) and the 

mass density p corresponds to a length scale in the Burgers equation. In this section, we recast the results of |^] in 
the language of ballistic aggregation of interest here 

We consider first the probability density of finding an interval [0, x) void of particles 



fioix) = / dyJ{y)J{x-y) 




( bx ( 1 \ I {w-\_-~W2)'^\ 



Ai(u;i)Ai(-u;2) 

which is plotted on Fig. ^. We have Yvaix^o l^o{x) — 1, and asymptotically for x ^ 00 



(102) 



/^o(x) = ^/ 7-^7- -j-^ 1 + O - ) ) • (103) 



[Ai'(-c.i)]' 




FIG. 4. Th e pr obability density ^io{^) to find no aggregates in an interval [0, x], for the parameters t = 1, p = 1 and P — 2 
(a = 1 in Eq.(Fo|)). 



The one point probability density ( |72| ) reads 



^In Burgers language, one speaks of the shock strength p/t and shock wavelength u — p/2 — ti]/p 

*In the present paper, the Brownian motion P(r) and the parabolic constraints /.^(r) have the sign opposite to the corre- 
sponding objects in namely —4>{y) and ~Su{y)- By invariance of Brownian motion under space reflexion, the functions / 
and J are the same in both papers, provided that M and P correspond to p and ~rj in the notation of Q. 
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P M 



^b'J T7-^ WJ ---- 



P M 



(104) 



2 J ^ ' \ M 2 

where X and J are given by Eqs.(b6Lfe8|). Integration over the momentum space yields the mass density distribution 



/ii(M) = 2b^MX{M)n{M) 



with 



whose asymptotic behavior reads 



dw- 



Hip) ^ 1, M ^ 0; ^(a*) ^ VTT&^^exp ^- 



~[2 



The shape of the mass density distribution (105) is plotted on Fig. ||. 



oo. 



(105) 



(106) 



(107) 




FIG. 5. The mass density fii{M) for the parameters t = 1, p = 1 and 13 = 2 {a — 1 in Eq.(10I:)) 



Its asymptotic behavior can be computed from (|100|) and (107) yielding the formulae 





and 



Ml 



(M) ~ 2V^b^/^M'/^ 



exp 



b^AP 
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uuibM 



M oo. 



(108) 



(109) 



The exact scaling function obtained here is quite different from a simple exponential exp(— M) suggested on the basis 
of numerical simulations in jlj. For example, one notices that small masses (M ^ t^/'^) are much more likely to be 
present in the system than suggested in ||l|, while large masses {M '3> t^^^) have a much smaller chance to be present. 
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Let us remark that the exponential form has been also found in [|j and ||] by solving the dynamical equations in a 
mean-field-like approximation scheme. In the framework of Burgers turbulence, more refined numeric al si mulations 
were performed by Kida where the small mass behavior of the scaling function compatible with Eq. f|l08|) has been 
found. Moreover, o ur r esult ( 108 ) is compatible with the rigorous upper and lower bounds derived in [|l2[| . The large 
mass behavior Eq.( |l09D fits into rigorous bounds of the type exp(— CM"^) found in and 
The density of particles of velocity V — P/M, given by 



dM MI{M)J(y - M/2)J{-V - M/2) 



(110) 



is plotted on FigJ 




FIG. 6. The velocity probability density ^i{V) for the parameters t — \, p — 1 and /3 = 2 (a = 1 in Eq.(llO)) 



From ( |110| ) and the asymptotic behavior ( [LOOD , (|101| ) of / and J, one can derive the bound 



(Ill) 



This shows that the large velocity behavior cannot be Gaussian and thus invalidates the numerically based hypothesis 
of Kida |]. 

The density of nearest neighbours (O) has the form 



/i2(l, 2) = 46^6'(X)5 {X + yf - 
&3 



X exp 



{{Y-f + (y+)3)) i{Ah)i{M2)j (-y+) J {Y-) 



(112) 



with X — X2 — Xi . Expression ( |112| ) permits to compute the collision frequency between the particles defined by 

Ai2(0,Pi,Afi,0+,P2,M2;t). (113) 



poo pOC 


Pi 


P2 


/ dPi / dP2 




' —00 'J —00 


Ml 


M2 



First, we notice that relation (|75|) implies the scaling 



V2{Mi,M2]t) 



iptr 
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j.2(M{,M^) 



(114) 



with 



/OO nOO 
-OO J — OO 



P]__P2_ 

Ml M2 



M2(l,2) (115) 



This last integral, upon inserting formula ( |112| ), gives 

i^2{Mi,M2) = 2bH4iM2{Mi + M2)I{Mi)I{M2)n{Mi + A'h) (116) 



with Ti given by Eq.(106). Clearly, the collision frequency does not factorize into a product of a function of Mi and 
a function of M2. This fact invalidates the weak mean- field hypothesis whose consequences where analysed in [|]. 

The two-point distribution function 

The scaling form p2(l,2) = pip2(2|l) of the two-point reduced number density of aggregates enables to study the 
long distance correlations in the system. The calculation requires the summation over all possible configurations of 
aggregated masses in between the particles 1 and 2 (formula ([43|) ) . In principle, this summation can be be performed 
as follows. Because of translation invariance, the conditional probability ^(2|1) 

IJl{X2,P2,M2\Xi,Pi,Mi) = {Pi,Mi\T{X2-Xi)\P2,M2) (117) 

can be considered as the kernel of an integral operator T{X) acting in the mass and momentum space. Then the 
series (^3|) can be summed up by applying the Laplace transformation. Setting 

1*00 

f{s) = / dXe'^''T{X) 
Jo 
1*00 

P2{s,PuMi,P2,M2)= / dXe-'''p2{Q,Pi,Mi,X,P2,M2) (118) 
Jo 

one finds by the convolution theorem 

00 

P2{S, Pi , Ml , P2, M2) = Pi {Pi ,Ml)J2{Pl: Ml |f '■(S)|P2, M2) 

= Pi(Pi,Mi)(Pi,Afi|f(s)(/-f(s))-i|P2,M2) (119) 

Inverting the operator / — T(s) and then performing the inverse Laplace transformation are not easy operations here. 
So, we shall proceed in a different way. In the context of the Burgers equation, ^2(1, 2) corresponds to the joint density 
to find two shocks at distance \X2 — Xi\. This density, in the Brownian interpretation, is the measure of paths that 
have contact points with two parabolas: one centered at Xi another at X2. Summing over all such paths amounts to 
sum over intermediate sequences of neighbouring aggregates in the particle language. We refer to |^| for details. One 
finds the following form of the two-point density 

P2{1, 2) = J {-Y+) /(Ml, Pi) [S {X + Y- - Y+) 

+9 {X + Y- - Y+) H {X, Y- , K+)] /(M2, P2) J (K,") (120) 

The 5— function term corresponds to the case when 1 and 2 are the nearest neighbours, whereas the function 
H {X, Yi, I2) embodies precisely the sum over all configurations of intermediate particles. The function H will not be 
reproduced here. The important fact is that its large distance asymptotic behaviour can be calculated leading to the 
folowing cluster property of the two-point function ||^ 

P2(Mi, Pi, M2, P2, X) - pi(Mi, Pi)pi(P2, M2) - -^'"/'^^ exp - bcjiX 



y 12 

X exp {-buJi{Yi- - y/)) J (-r+) I{Mi)I{M2)J {Y2) , X -> 00. (121) 

We see that although the particles in the initial state are not correlated and the motion between collisions is free, the 
aggregation process induces dynamic correlations in the course of time. However, these correlations have a very short 
range since they stay dominated by the rapidly decaying cubic exponential factor exp(— 6''X''/12). 
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X. RIGOROUS SOLUTION OF THE DYNAMIC HIERARCHY 



We have shown in Section IV how the infinite hierarchy (|T]) describing the aggregation process could be reduced 
to a system of two coupled equations (^3|), (^J), satisfied by the one-particle density fj,i{l;t), and the conditional 
probability density fj,(2\l;t). It is the right moment now to show that the constructed densities ([72|), ( |73| ) provide a 
rigorous solution to the aggregation dynamics. 

In ( [72I) and (73) the densities and ii{2\l;t) are expressed in terms of the probability weights I{M,P) and 

J{Y) (to recover the time variable t the scaling relations (^^ have to be used). In order to evaluate the rate of 
change of the state of the aggregating gas it is thus sufficient to calculate the time derivative of functions I{M, P; t), 
J{-Y+;t) and J(r";tV where = tP/M ± M/2p. A straightforward calculation can be performed starting from 
the defining formulae (|5^), (js^), whose explicit form can be found in |^]. Technically it is simple but lengthy, so we 
shall not reproduce it here. Taking then the continuum limit one finds the following system of equations 



d_ 



I{M, P\t)= / dMi / dPi / dAh / dP^ 



xS 



Ml 



Po 







M' 


\Mi 


M2) 





M2 

I{Mi,Pi-t)I{M2,P2;t) 



5{M - Ml - M2)S{P -Pi- P2) 



(122) 



d_ ( P_ M 

dt \M^~Yp'\ 



dMi 



X , Pi Ml 
xIiPi,Mi;t)j[-t--;t 



P 
M 



Pi_ 

Ml 



M + Mi 



2p 



(123) 



A somehow laborious but straightforward analysis permits to check that the above system of equations, when 
combined with the relations ([7^), (|7^), implies the fundament al dy namic equations (p^), (|3^). As one could expect 
from the interpretation of functions / and J given in Section VII, the time derivative oiI(M,P;t) is positive and 
produces gain terms in the kinetic equation (33), whereas the negative time derivatives of factors J{Y^; t) and J{Y~; t) 
generate loss terms. 



XI. CONCLUDING COMMENTS 

The hierarchy equations ( |2l| ) describe the aggregation dynamics in a one-dimensional gas for arbitrary initial 
conditions. We considered here in detail the case where statistical correlations existed only between the states of 
the nearest neighbours. Under the precise condition (|3^), the evolution preserved this property, and the state of the 
system remained entirely characterized by the number density of aggregates and by the conditional distribution of 
nearest neighbours. This permitted to reduce the hierarchy to two coupled equations. 

Rather than solving the obtained system of equations we determined the relevant distributions on the basis of their 
definitions, and only a posteriori we could verify that they satisfied the dynamical laws. So, finding a systematic way 
of solving the system (p3|), ( ^^ remains an open problem. 

When discussing the continuum limit we started from a very simple discrete initial condition, where no correlations 
were present at all. In the course of time, the system turned out to build up correlations by its internal dynamics, 
but only between the pairs of nearest neighbours. It is thus tempting to conjecture that the factorized form ( p2[ ) of 
the many particle densities represents the asymptotic long time structure of the system for general initial conditions, 
as only the nearest neighbours get corr elate d by the process of aggregation. Let us recall that the correlations are 



rapidly decaying with the distance (see ( 121 ) 



A comment concerning the continuum limit seems also quite important. It has been noted at the end of Section VII 
that introduction of properly scaled variables permitted to interpret the continuum limit as a long time limit. One 
can thus expect that the dynamical properties derived here are quite universal, and do not depend on our particularly 
simple choice of the initial distribution. In particular, it would be interesting to clarify how large is the class of initial 
states which leads for long times to the mass density distribution plotted in Fig.(||). 

Let us finally notice, that the dynamical scaling of the mass distribution predicted on the basis of intuitive arguments 
in 1^ agrees with the rigorous scaling relations (]59|). The prediction of the Brownian motion exponent a = 2/3 in |^ 
seems independent of the initial condition, whereas we derived it starting from a particular state. It would be thus 
interesting to clarify to what extent this type of scaling is universal. Indeed, our analysis showed that the aggregation 
dynamics was in principle compatible with other values of a, not necessarily equal to 2/3. 
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APPENDIX A: 



The moments of the distribution Gk (pq) can be obtained from its generating function 7(s,^) 



fc=l 



€=0 



with 



k=l -^0 



(Al) 



(A2) 



This generating function can be computed by applying the Sparre- Andersen theorem to the process with independent 
increments governed by the Gaussian distribution (tt)^^/^ exp(— g^) (section XVIII. 3 of jl^) 



7(s, ^) = exp 



„fc /-oo 



.k=l 



(A3) 



Using ( |A3| ) we find 



^.'^ (1),= 7(^,0) = exp ^ 



oo 



and thus 



(1), = 



For the first and second moments we find 

oo 
k=l 

and 



r(fc + i/2) 
V^r{k + 1) 



\k=l 



'irk 



2k / 



, k — > oo. 



.^0 " ^1 2V^ 



(A4) 



(A5) 



(A6) 



k=l 



de 



5=0 4yr^ 



oo t. 



(A7) 



To proceed, we need the following Tauberian theorem (Theorem 5, section XIII. 5 in |Tj]): 

Let Cfe > be a monotonic sequence and suppose that the series X^feli ^feS*^ converges for < s < 1. Then, if L 
varies slowly at infinity ^ and a > 0, each of the two relations 



E 

k=l 



CkS 



' -l' ' 



(1-.) 



1 - s 



s 1 



(A8) 



varies slowly at infinity if limt-,oo L{tx) / L{t) = 1 for all x, a condition which is verified in the subsequent applications 
where L is a constant. 
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implies t he o ther. 
From ([Aq), we have 



and thus 



For the second moment we find 



leading to 



Ck- T^^k'^-'Lik), fc->oo (A9) 
T(a) 



OO ^ 

fe=i ^ ' 



k^^. (All) 



OO 



fc=i 



/c->oo. (A13) 



To control the remainder in ( |90D by a limited Taylor expansion up to second order in q we use the fact that Gk (q) 

is non negative and that the second g-derivative of e~'^^/^^'^ K,n,iy{MQ, vMq + y^^g, M, P) has a g-integrable bound 

for Af > uniform when m and Afo are in a neighbourhood of zero (irrespective of the sign of v) . This is indeed the 
case because one can check that Ku{Mq, Pq, P) as well as its first and second Po-derivative obey Gaussian bounds 
of the form Ci exp(— C2P0 ) with Ci and C2 independent of m and Afo, < M 
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